A Phenomenological Theory of The Pseudogap State 
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An ansatz is proposed for the coherent part of the single particle Green's function in a doped 
resonant valence bond (RVB) state by, analogy with the form derived by Konik and coworkers for 
a doped spin liquid formed by an array of 2-leg Hubbard ladders near half-filling. The parameters 
of the RVB state are taken from the renormalized mean held theory of Zhang and coworkers for 
underdoped cuprates. The ansatz shows good agreement with recent angle resolved photoemission 
(ARPES) on underdoped cuprates and resolves an apparent disagreement with the Luttinger Sum 
Rule. The transition in the normal state from a doped RVB spin liquid to a standard Landau 
Fermi liquid, that occurs in the renormalized mean held theory, appears as a quantum critical point 
characterized by a change in the analytic form of the Green's function. A d-wave superconducting 
dome surrounding this quantum critical point is introduced phenomenologically. Results are also 
presented for the Drude weight and tunneling density of states as functions of the hole density. 

PACS numbers: 74.20.Mn, 74.25.Jb, 79.60.-i 



I. INTRODUCTION 

The cuprate superconductors have attracted enormous 
interest not just because of their high temperature un- 
conventional superconductivity but also because of their 
highly anomalous properties in the normal phase (for a 
review see T. Timusk and B. StattA). These show strong 
deviations from standard Landau Fermi liquid behavior 
particularly in the underdoped region where the devia- 
tions are most evident in the pseudogap (or spin gap) 
phase. Among the most spectacular of these deviations 
is the observation of a Fermi surface in photoemission 
(ARPES) experiments that consists, not of a closed con- 
tour, but only of 4 disconnected arcs centered on the 
Brillouin zone diagonals. 2 ^ Pseudogaps in the single par- 
ticle spectrum truncate the Fermi surface at the saddle 
points. In addition the density of charge carriers is deter- 
mined not by the conduction electron density but by the 
hole density doped into the stoichiometric parent Mott 
insulator A These facts indicate that the pseudogap phase 
should be viewed as an anomalous precursor to the stoi- 
chiometric Mott insulator. 

Very soon after the discovery of the high temperature 
superconductors and before these anomalous pseudogap 
properties were measured, Anderson proposed that the 
key to understanding these unique materials lies in what 
he called Resonant Valence Bond (RVB) behavior^ He 
proposed a description based on lightly hole doped spin 
liquid of spin singlets. Rather than forming a hxed ar- 
ray of singlets, strong quantum fluctuations among the 
antiferromagnetically coupled S = 1/2 spins lead to a 
superposition of singlet configurations, i.e. the bond sin- 
glets resonate between many configurations. This ele- 
gant concept explains many key features of the pseudo- 
gap phase as emphasized in the recent review by Ander- 
son and coworkers^ But the strong onsite correlations 
among the electrons in this precursor to a Mott insulator 
are difficult to treat analytically. This has hindered the 



development of complete theory for the RVB phase and 
the cuprates in general. 

In this recent review^, Anderson and coworkers point 
out that the early renormalized mean field theory (RMF) 
introduced by Zhang and coworkers 6 for the RVB 
phase predicted many of its key features using simple 
Gutzwiller renormalization factors to describe the strong 
correlations. In the intervening years much progress has 
been made on a first principle treatment of these correla- 
tions based on gauge theories (for a review see Lee, Na- 
gaosa and Wen—) . The functional renormalization group 
(RG) approach developed by Honcrkamp and coworkers 
has also given new insights starting from weak coupling. 8 
However, we still lack even a consistent phenomenological 
description of this anomalous pseudogap phase that ties 
together key features analogous to the Landau theory for 
standard Fermi liquids. The purpose of this paper is to 
take a step in this direction. To this end we introduce an 
ansatz for the form of the self-energy and thereby for the 
form of the single particle Green's function. In this we 
are guided by the recent theory by Konik, Rice and Tsve- 
lik (KRT) for the form of the Green's function in a doped 
spin liquid consisting of an array of 2-leg Hubbard lad- 
ders coupled by long range inter-ladder hopping. 9 By a 
careful choice of the form of this hopping, KRT could jus- 
tify a random phase approximation (RPA) to obtain the 
two dimensional Green's function from the single ladder 
form. They showed that for a light doping away from 
one electron/site a novel form of the Green's function, 
G(k,u>) resulted. In particular they showed that the be- 
havior of G(k, 0) which enters the Luttinger Sum Rule 
(LSR)—, was quite distinct from that of a normal Lan- 
dau Fermi liquid. The LSR relates the area enclosed by 
the contours where G(fc,0) changes sign from positive 
to negative to the total electron density. In a Landau 
Fermi liquid, sign changes in G(fc,0) occur only through 
an infinity in G(k, 0) at a closed Fermi surface. In the 
doped spin liquid, KRT found that the sign changes oc- 
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cur not only through infinities on Fermi pockets but also 
through zeroes lying on a separate surface which they 
call the Luttinger surface of zeroes. The Fermi pockets 
enclose an area determined by the hole density while the 
Luttinger surface encloses an area given by the density of 
one electron per site in the parent stoichiometric ladders. 
This behavior followed from the form of the self energy 
in the coherent part of the single ladder Green's func- 
tion and its modification through the inter-ladder hop- 
ping processes. The single ladder Green's function was 
derived by Essler and Konik only in the weak coupling 
limit^i However, the fact that weak and strong coupling 
are continuously connected in ladders and that the form 
is qualitatively similar to that obtained numerically by 
Troyer et al in strong coupling t — J ladders 12 , leads us to 
believe that this general form is characteristic of doped 
spin liquid and is not restricted to weak coupling. 

In section 2 we introduce our ansatz for the form of the 
self energy for a two dimensional lightly doped RVB spin 
liquid. Our ansatz is based on a generalization of the 
KRT form for the doped spin liquid formed in an array 
of 2-leg Hubbard ladders^ It also contains a Luttinger 
surface of zeroes which encloses a commensurate density 
in addition to hole Fermi surface pockets. In this case 
the Luttinger surface coincides with the so-called umk- 
lapp surface introduced by Honerkamp and coworkers 8 
as the surface where umklapp scattering processes ap- 
pear to open up a charge gap in the weak coupling RG 
flow. 

It is also interesting to note that this Green's function 
is closely related to the form recently proposed by Tai- 
Kai Ngi£. He started from the strong coupling limit and 
introduced spin-charge separation by factorizing single 
electron operators into a product of spinon and holon op- 
erators. Following earlier work by Wen and Leei^ he in- 
troduced a phenomcnological attraction between spinon 
and holon which leads to a binding between spinon and 
holon. He obtained a form for the Green's function which 
also contains a coherent quasiparticle pole and has close 
similarities to our ansatz. 

In section 3, we analyze the consequences of our ansatz 
for a variety of properties. First we obtain the hole den- 
sity dependence of our psuedogap and other key parame- 
ters from the RMF results of Zhang et al. 6 Based on these 
choices we obtain estimates of the hole density depen- 
dence of many observables, e.g. the Fermi velocity vf and 
wavevector kp along the Brillouin zone (nodal) direction, 
the minimum gap contour near the saddle point (antin- 
odal direction), the tunneling density of state (DOS), the 
Drude weight, and the shape and form of the hole Fermi 
pockets. From the phenomenological form of the Green's 
function we obtain the quasiparticle dispersion and spec- 
tral weights that characterize the coherent part of the 
single electron Green's function. A key point in our phe- 
nomenology is that there is a critical hole concentration 
above which the spin liquid anomalous self-energy van- 
ishes. As consequence there is a form of quantum critical 
point (QCP) which separates two topologically distinct 



forms for the Green's function. Below the critical hole 
density, G{k, 0) is characterized by coexisting Luttinger 
surfaces of zeroes and hole pocket Fermi surface of infini- 
ties while above the critical concentration G(k,0) dis- 
plays only a closed Fermi surface of infinities as usual in 
a Landau theory. Many features of this QCP resemble 
those deduced by Loram and coworkers from an analysis 
of a variety of experiments^. 

In section 4, we examine how this Green's function is 
modified when the system enters a d-wave superconduct- 
ing state. This again is introduced phenomenologically 
and no attempt is made here to derive the parameters 
of the superconductivity although we do offer a number 
of reasons why the hole pockets should have a d-wave 
superconducting instability. 

Lastly section 5 is devoted to further discussions and 
conclusions. 



II. ANSATZ FOR THE SELF ENERGY IN A 
DOPED RVB SPIN LIQUID 

A doped RVB spin liquid has many properties which 
are quite distinct from the standard Landau Fermi liquid 
which follows from treating the interactions in perturba- 
tion theory. Moreover these anomalous properties cannot 
be ascribed to a broken symmetry or to the appearance 
of a new order parameter. As a result it is a challenge to 
construct a consistent theory even on the phenomenologi- 
cal level of a doped RVB spin liquid. We start with a brief 
review of the recent work by KRT— who could obtain the 
form of the single particle Green's function G(k,ui) in a 
doped spin liquid. In particular, they examined the form 
of the LSR which applies to the zero frequency Green's 
function G(fc, 0) in the doped spin liquid. The key point 
about the LSR as emphasized in the famous textbook 
by Abrikosov, Gorkov and Dzyaloshinskii (AGD)i& and 
more recently by Tsvelikii, is that derivation of the LSR 
is very general and is not limited to perturbation theory. 
The LSR relates the total electron density, p, to the area 
in the fc-space where G(k, 0) > 0. In two dimensions it 
takes the form 

P=lha\ d " k (!) 
( 27r r Ja(k,o)>o 

An important point that these authors emphasized is 
that the sign change from positive to negative values of 
G(k,0) is not restricted to an infinity in G(k,0) such as 
occurs at the Fermi surface of a Landau Fermi liquid. It 
can also occur through a zero in G{k, 0) as, for example, 
in the case in the BCS theory of superconductivity. KRT 
considered a doped spin liquid consisting of an array of 
two-leg Hubbard ladders. At half-filling in a single ladder 
the Fermi surface consists of four points without interac- 
tions, but it is completely truncated when the repulsive 
interactions which lead to both spin and charge gaps, 
are introduced. All spin and charge correlation functions 
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are strictly short range so that this system is a true spin 
liquid. 

An explicit form for the single particle Green's function 
has been derived by Konik and coworkers in the limit of 
weak repulsion^i*i& Around each of the four Fermi points 
G(k,u>) can be split into coherent and incoherent parts, 



Zgiu + £q(fc)) 

lo 2 - el(k) - A 2 



(2) 



Here e a (k) is the bare dispersion near the corresponding 
Fermi wavevector kp a , z a is quasiparticle weight ~ 1, 
and A the single particle gap. Note while the coherent 
part has a form similar to the diagonal Green's function 
in BCS theory, there is no off-diagonal component of the 
Green's function in this case. This general form is also 
compatible with the numerical results of Troyer et al for 
the strong coupling limiti^. G^(k,0) changes sign from 
positive to negative through a zero in G^(k, 0) which oc- 
curs when e a (k) passes through zero as k crosses kp a , 
and the area with G^(k, 0) > is unchanged by the in- 
teractions, satisfying the LSR. 

Starting from this result, KRT could derive G^(k,Lo) 
for an array of 2-leg Hubbard ladders coupled by a par- 
ticular form of long range inter-ladder hopping chosen so 
that it could be treated using a random phase approxi- 
mation (RPA) 



G 



RPA 



(k,Lu) = l/{G^k,Luy 1 -t ± (k ± )} (3) 



where k — [k, kx.)- They showed that this form, for 
values of t± large enough, leads to the appearance of 
electron and hole pockets distinct from the lines of ze- 
roes in G RPA (k,0). The zero contours remain as lines 
at k = kp a , independent of the transverse component of 
the wavevector, kj_- In the presence of hole doping, the 
hole pocket expands and the electron pocket shrinks. In 
this doped spin liquid G B - PA (k, 0) has sign changes both 
through infinities along the Fermi pockets and through 
zeroes located on the lines k — kp a . Note although 
the zeroes lines appear at incommensurate wavevectors 
in general and not at the Brillouin zone as in a stan- 
dard Bloch insulator, the total area they enclose is com- 
mensurate and equals 1 electron / site, irrespective of the 
inter-ladder hopping strength. The final form obtained 
by KRT for the coherent part is 



G 



RPA 



lu - e a {k) - t±{kx) - A 2 /(^ + e a (k)) 



(4) 

This can be interpreted as a ladder self energy Si(fc, lu) — 
A 2 / (ui+e a (k)), where e a (k) — at the fc-points where the 
gap opens up in the parent insulating two-leg Hubbard 
ladder array. 

In this work we consider doping a two-dimensional res- 
onant valence bond insulator. We start from the RMF for 
such an insulator. This approximation treats the effect of 
strong correlations through renormalization factors cal- 
culated by a description that goes back to early work by 



Gutzwiller— . The renormalized t — J Hamiltonian takes 
the form of an effective Hamiltonian 6 



H, 



eff 



9tT ■ 



g s J ^2 S i 



(5) 



with the kinetic, T, and spin energy terms modified by 
factors gt and g s respectively 



9i 
9s 



2x 
1 + x' 
4 

(l+.x) 2 



(6) 



for a hole doping of x per site. At half-filling g t = 
leaving only the magnetic energy. The RVB ansatz 
factorizes the spin energy introducing both Fock ex- 
change, Xi,j = ( c l,<r c j>) and pairing, A itj = (cj >CT Cj jCT ) 
expectation values. The factorization procedure is not 
unique but the spin quasiparticle dispersion that results 
is unique, = (3g s J/8)(cos 2 fc^+cos 2 ky) 1 / 2 . Upon hole 
doping gt > 0, and coherent quasiparticle poles with a 
small weight g t appear in the single particle Green's func- 
tion. By analogy with the KRT form for the doped spin 
liquid discussed above, we make the following ansatz for 
the coherent part of G(k, lo) in a doped RVB spin liquid 



G 



RVB 



(*,«) 



9i 



u-m-K/^ + Uk)) 



Gi 



(7) 



where k = (k x , k y ), 



£ (fc) = — 2t(x)(cosk x + cos k y ) 
Afl(fc) = Ao(x)(cos k x — cos k y ) 
£(fe) = £o(k) — 4t'(x) cos k x cos k y 

-2t"(x)(cos2k x + cos2k y ) - fi p (8) 

Eq0 is analogy to EqQ] for the coupled ladder system, 
and £(fc) — £, (k) is analogy to t±(k±) in EqQ] In the 
renormalized dispersion we include hopping terms out to 
3 rd nearest neighbor with coefficients 



t(x) 

t'(x) 
t"{x) 



9t(x)t 

9t(x)t' Q , 
9t(x)t'o 



-g s (x)Jx, 



(9) 



The RVB gap magnitude function Ao(x) is also taken 
from the RMF theorji 6 .. The parameter \i represents a 
shift of the energy band so that the chemical potential is 
always the zero of the energy. We determine /j, p from the 
LSR on the total electron density. 

First we consider the limit of zero doping, x — > 0. 
In this case gt(x) — > and the quasiparticle disper- 
sion reduces to the spinon dispersion and the quasipar- 
ticles have the vanishing weight in this limit in the sin- 
gle particle Green's function, G(k,oj). At small but fi- 
nite x the zero frequency Green's function G RVB (k,0) 
that enters the LSR has lines of zeroes when ^,o(k)(— 
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— 2t(x)(cos k x + cosky)) = 0. The Luttinger contour of 
zeroes in G RVB (k,0) then consists of straight lines con- 
necting the points (±7r,0) and (0, ±tt). This Luttinger 
contour coincides with the antiferromagnetic (AF) Bril- 
louin zone. But more relevantly it coincides with the 
umklapp surface which appears in functional RG calcu- 
lations on the weak coupling 2D t — t' — U Hubbard 
model. Honerkamp8>, Laeuchli and coworkers 2 ^ found 
that umklapp scattering processes in both particle-hole 
and particle-particle channel which connect points on this 
surface grew strongly at low energies and temperatures 
leading them to propose that an energy gap would open 
up on this surface below a critical scale. Further because 
of a close analogy to the behavior found in the case of 
a half-filled 2-leg Hubbard ladder, they proposed that 
this gap was not driven by long range order, but rather 
was a sign that a RVB spin liquid with short range order 
would form below a critical scale. Thus our ansatz for the 
Green's function is fully consistent with these proposals. 

A second feature that follows from our ansatz for 
G RVB (k,w) is the appearance of hole pockets at finite 
hole doping. These will be reviewed in detail in the next 
section. The hole pockets define Fermi surfaces where 
G RVB (k,0) changes sign through infinities and contain 
a total area equal the hole density. The LSR is satis- 
fied since the area with G RVB (fc,0) > is bounded by 
the Luttinger surface which contains 1 el/site, minus four 
hole pockets which have a total area related to the hole 
density. 

The phenomenological form G RVB (k,0) for a hole 
doped RVB spin liquid can be straightforwardly general- 
ized to a d-wave superconducting state 

= <fc/[w-£(fe)-E fl (fc,w) 

- |A s (fc)| 2 /( w + £(fc) + £ fl (fc, -«))] (10) 

where Sfl(fe, w) is the RVB spin liquid self energy from 
Eq.0 

^(fc,c) = |A fl (fc)i 2 /(^+afc)) (ii) 

and As(fc) is the d-wave superconducting gap function. 
To analyze the LSR on Gf o/l (fc,0) we note first that 
Sfl(fe,0) — > oo on the surface where £ (fe) = 0. Thus 
Gf o/l (fc,0) continues to have a Luttinger surface of ze- 
roes on the umklapp surface as in the normal RVB spin 
liquid. However, there is now an additional set of Lut- 
tinger surfaces defined by the contour that satisfy 

£(fc) + £ fl (fc,0) = 

But these are just the Fermi surface of the four hole pock- 
ets in the normal phase which have now been converted 
to a Luttinger surface of zeroes in the superconducting 
state. Thus the form Eq. ^U] continues to satisfy the LSR. 

Lastly we remark that along the Brillouin zone diag- 
onals both Afi(k) and Ag(fc) vanish in our phenomeno- 
logical form. As a result exactly along these directions 



there is only a single quasiparticle pole which crosses 
the Fermi energy at a Fermi wavevector determined by 

£(k)\k=(k F ,k F ) = 0- 

Our phenomenological form for G RVB in the normal 
and superconducting phases are similar but not identi- 
cal to several other recent proposals. It is very close to 
the form Tai-Kai Ngi^ derived based on spin-charge sep- 
aration but with an added phenomenological attraction 
between spinon and holon which leads to binding and 
therefore to quasiparticle poles and a self energy form 
similar to Eq.^l]. He, however, restricted his analy- 
sis to the case of only nearest neighbor hopping. Ear- 
lier Norman and coworkers from an analysis of ARPES 
data on BSCCO samples around optimal doping intro- 
duced a d-wave superconducting self energy modified to 
included both normal state single particle scattering and 
pair scattering^! Recently Honerkamp2£ has speculated 
on an adaptation of the form introduced by Norman et 
al to describe a truncation of the Fermi surface through 
the opening of energy gap in the antinodal regions near 
(±7r, 0), (0, ±7r) as a phenomenological description of the 
RVB spin liquid. 



III. ELECTRONIC PROPERTIES OF THE 
NORMAL PSEUDOGAP PHASE 

In this section we discuss the electronic properties that 
follow from our phenomenological form for the Green's 
function and compare these to experiments on the normal 
pseudogap state. 

We begin with the LSR shown in Eq.pQ. In Fig^a- 
e) we show that the contours on which G RVB (k,u> — 
0), defined in Eq.[7|, changes sign. We chose val- 
ues for hopping parameters in Eq.[5], appropriate to 
Ca2- x Na x Cu02Cl2- The underlying band structure pa- 
rameters were obtained by a tight-binding fit to the anti- 
bonding 3d x 2_ y 2 —2p x ( v j band calculated by Matthciss 23 . 
These were then renormalized by the Gutzwiller factors 
defined in Eq.jfj], leading to the values shown in FigEl 
The sign changes in G RVB {k,u> = 0) occur on a Lut- 
tinger surface of zeroes, which coincides with the umk- 
lapp surface, and on a Fermi pocket of infinities. The 
parameter fi p was adjusted at each x to give the correct 
area for the hole pockets (see FigEJ - 

As we can see from Fig^ the hole pocket evolves 
gradually into a more normal surface in panel(e) as x 
increase. The spectral weight of the quasiparticle pole 
varies strongly around the pocket as illustrated in Fig|3] 
In particular it is very small on the outer edge of the 
pocket closest to the Luttinger surface of zeroes and van- 
ishes completely along the nodal directions ((1,1) direc- 
tions). Along this direction there is only a single sign 
change in G RVB (k, to = 0) at the inner pocket edge. 

In FigQ] the dispersion of the quasiparticle poles Ei(k) 
in the coherent part of G RVB together with their spec- 
tral weight, Zi(k), are shown. We see that in general 
there are two quasiparticle bands Ei{k) with strongly 
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(a): x=0.05 



(b): x=0.10 



(c): x=0.14 




FIG. 1: (Color online) Contours on which G(k,0) changes 
sign at various hole concentrations x are shown in (a)-(e). In 
the shaded area G(fc,0) > 0, satisfying the Luttinger Sum 
Rule. In the normal pseudogap phase, the line connecting 
(tt, 0) — (0, 71") is the Luttinger surface of zeroes and the pockets 
in the thick line represent the infinities of G(k, 0). The values 
of the parameters used here are given in Fig|5| The evolution 
of the contours of infinities in G(k, 0) is illustrated in (f). 
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FIG. 2: (Color online) The values of the parameters t(x), 
t'(x), t"(x), n p and A(x) that enter in Eq. 18191 used in 
the present calculations. (\ = 0.338, J /to = 1/3) Re- 
sults presented in the text are for a choice of t' /to = —0.3, 
to /to — 0.2, estimated from the calculated band structure of 
Ca 2 Cu0 2 Ch& 
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(b):x=0.10 
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(e): x=0.20 
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FIG. 3: (Color online) The spectral weight distribution 
z(k)/gt around the hole quasiparticle pockets in the normal 
pseudogap phase. 



varying spectral weight. Over most of the Brillouin zone 
the bands are separated by an energy gap with only the 
lower band occupied. The distribution of spectral weight 
between the two bands is determined by the proximity 
to the umklapp surface. Near this surface the weight is 
equally divided but away from this surface only a small 
admixture is induced by the anomalous self energy 
Also along the nodal directions — > 0, so that in these 
directions the two bands coalesce into a single band. 

In order to compare to ARPES experiments on 
Ca2~ x Na x Cu02Cl^, we prepared a single plot which 
combines contours of the hole Fermi pockets and also the 
minimum energy gap lines in other parts of zone. These 
are illustrated in Fig[S]together with the ARPES results. 
There is good agreement between the two sets of curves 
with the exception of the outer edge of the hole pockets 
which has not been reported in the ARPES experiments. 
The predicted spectral weight of the quasiparticle band 
on these outer edges is very small (see Fig|3J). Nonethe- 
less it would be important to search more closely for a 
weak signal on this edge. In FigHjl we plot the dispersion 
of the lower quasiparticle band E(k) on the "Fermi sur- 
face" and on the hole pockets for several dopings. The 
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FIG. 4: (Color online) Quasiparticle dispersion in the normal 
pseudogap phase (Eq.Q) along some high symmetry direc- 
tions. The thickness of the lines is proportional to the spectral 
weight z(k)/gt of the quasiparticle. 



angular dependence of |i?(fe)| on the "Fermi surface" is 
shown in panel (e). The dispersion drops to zero rather 
sharply as the angle touches the pocket direction. 

In Fig[7| we show the hole density dependence of the 
two features, the Fermi wavevector, kp (with kp = 
fcp(l, 1)) at the hole pocket along the nodal line (1,1) and 
the intercept of the minimum energy gap line at the Bril- 
louin zone boundary connecting (tt, 0) — {it, 7r) antinodal- 
kp (with kp = (kp,Tr)). Again the agreement is very 
good with ARPES experiments on Ca^-x^ a-xCuO^Cl^. 
These two wavevectors have been interpreted as key pa- 
rameters of an underlying Fermi surface. But such an 
interpretation clearly violates the LSR since the enclosed 
area would require electron doping in contrast to our 
ansatz for the Green's function. As we remarked ear- 
lier, our ansatz reconciles the ARPES measurements and 
the LSR. 

Also in Fig[7|we show the Fermi velocity vp along the 
nodal direction. This shows rather more variation with 
the hole doping, x, than the direct calculations from the 
Gutzwiller projected variational wavefunctionSi but the 
qualitative behavior is similar. However, it deviates sub- 
stantially from experiment with increasing x. The dis- 
crepancy of vp at large x could be due to the oversimpli- 
fied model. 

Another quantity of interest is the coherent quasipar- 
ticle contribution to the tunneling density of state (DOS) 
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FIG. 5: (Color online) Comparison between our theory and 
some recent ARPES experiments on Ca2- X N ' aCuO-zGl by K. 
M. Shen et al^ The experimental results are re-plotted in 
panels (A-F). Panels (A-C) show the distribution of spectral 
weight in the Brillouin zone within a ±10meV window around 
the Fermi level. The open/solid circles in panels (D-F) are de- 
tected by this experiment to determine the Fermi surface. The 
pockets in panel (A-I) show the infinities of G(k, 0) in the nor- 
mal pseudogap phase. The "Fermi surfaces" (blue and black 
curves) shown in panels (D-I) are defined by the minimum dis- 
tance from the lower quasiparticle band to the Fermi level in 
our calculation along radial directions centered at (n, n). The 
thickness of the curve in panels (D-I) represents the spectral 
weight z(k)/g t of the quasiparticle. 



defined as 



N t (oj) 



J2 Zi{k)5((j - Ei(k)) (12) 



fc,i=l,2 



The value at the chemical potential, Nt(0) is shown 
for a series of hole densities in Fig|Sfa). iVx(O) rises lin- 
early as the hole density x is increased from zero in the 
pseudogap phase. This behavior is similar to the Drude 
weight, D^, calculated by integrating around the Fermi 
surfaces of the hole pockets 



D ' lv ~ (2tt) 2 h 



dSf 



(13) 



In the square lattice D^ v is a diagonal tensor with ele- 
ment D^^ shown in FigOJb). 

The density dependence of the RVB gap parameter, 
Ao(x) was taken from the RMF results of Zhang et al£ 
This gives a linear drop in Ao with Ao(x) — > 0, as x — > 
x c (see in FigJ2{e)). We have chosen a value of x c = 
0.2. This linear drop in Aq(x) as x increases is in line 



7 
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(b): x=0.10 
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FIG. 6: (Color online) Panel(a-d) show the dispersion of the 
lower quasiparticle band E(k) on the "Fermi surface" and on 
the hole pockets for various dopings. In panel(e) we show the 
angular dependence of \E(k)\ on the "Fermi surface". (j> is 
defined in the inset of panel (e). 
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FIG. 7: (Color online) Nodal Fermi wavevector Jcf, 'antin- 
odal Fermi wavevector' hp, Fermi velocity vf (with values 
to = 0.5eV, lattice constant a = 4. OA). The experimental 
data in panel (a,b) for nodal and 'antinodal Fermi wavevector' 
are for Ca2- X N a x Cu02Cl^ while the experimental results in 
panel (c) are for Ca2- x Na x Cu02Ch, La2- x Sr x CuC>4, and 
Bi2Sr 2 CaCu20 fi+ ^. 



with ARPES experiments as emphasized in the recent 
review by Anderson et al£ The vanishing of Ay (a?) at 
x = x c causes a change in the form of G(k,u>) with a 
transition from the doped spin liquid RVB character to 
that of a standard Landau Fermi liquid. Then contours 
of sign changes in G(k, 0) that enter the LSR change their 
topology and the Luttinger surface of zeroes disappears. 
Only a standard closed Fermi surface of infinities exists 
at x ^> Xc • 

The detailed form of the transition from a doped RVB 
spin liquid to a Landau Fermi liquid depends on the band 
parameters. For the choice that we made to correspond 
to Ca-i-xN a x CuOiClii there is an additional topologi- 
cal change in the LSR contours clS X — ^ Xr . As illustrated 
in Fig^ for a narrow range of x < x c a new set of elec- 
tron Fermi pockets appear close to the saddle points, out- 
side the Luttinger surface of zeroes. These new electron 
pockets merge with the hole pockets inside the Luttinger 
surface of zeroes at x — x c to give a Fermi surface that 
crosses the umklapp surface at x > x c . In our parameter 
choice the saddle points in the band structure are occu- 
pied for a range of hole densities for x > x c which requires 
that the Fermi surface cross the umklapp surface in this 
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FIG. 8: (Color online) Density of states at the Fermi level 
Nt(0) and the Drude weight D M/J in the normal pseudogap 
phase. 
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FIG. 10: (Color online) Regular and anomalous Green's 
functions G(k,uj), F^(k,uj) which appear in the coupled 
En. |lillS| . 



tively similar to that inferred by Loram, Tallon and col- 
laborator from experiments^***— Their analysis empha- 
sized the singularities in thermodynamic quantities such 
as specific heat, magnetic susceptibility etc. Our phe- 
nomenological theory is restricted at present to zero tem- 
perature and so is not suitable for detailed comparison. 
However, we note that the main conclusion they drew, 
that a partial gap opened up at x < x c , is in agreement 
with our phcnomcnological theory 



FIG. 9: (Color online) Density of states Nt(cj) in the normal 
pseudogap phase. The zero energy values Nt(0) are shown in 
FigUa). 



density range. For this parameter choice, there are two 
closely spaced topological changes in the LSR contours. 
The first is a Lifschitz-type transition in which the shape 
of the Fermi surface of infinities undergoes a topological 
change. The second at quantum critical point 

(QCP) associated with opening of a single particle gap at 
the onset of the RVB spin liquid for x < x c leading to a 
Luttinger surface of zeroes coinciding with the umklapp 
surface. 

These two transitions cause singularities in the density 
dependence of the DOS and the Drude weight as illus- 
trated in FigJSJ The strongest singularity appears at the 
QCP at x = x c - A substantial jump appears in the DOS 
at x = x c . The DOS continues to rise as x increases 
beyond x c which is associated with the approach of the 
Fermi energy to the van Hove singularity at the saddle 
points in the band structure. A divergence of the DOS 
at the Fermi energy at x > x c has not been observed 
to our knowledge but this might be due to different band 
parameters or possibly to a suppression of the divergence 
due to impurity scattering. The full energy dependence 
of the DOS is shown in Fig|5] 

The QCP in our phenomenological theory is qualita- 



IV. ELECTRONIC PROPERTIES OF THE 
SUPERCONDUCTING STATE 

We turn now to the evolution of the electronic prop- 
erties when the system enters a d-wave superconducting 
state. Superconductivity is introduced by the addition 
of the new term in the self-energy which follows from the 
standard Green's function theory of a superconductor—. 
Solving the coupled equations which connect the regular 
(G s (k, u>)) and anomalous (F(k,u>)) Green's functions 
(see FigllDJ 

(co-£(fc)-£ K (fc,^))G 5 (fc, W ) 
-»A s (fc)F t (fc,w) = l (14) 

(-w-e(fc)-EH(fc,-«))Ft(fc sW ) 
-iAl(k)G s (k,oj) = (15) 

leads the result quoted earlier in Ea.|l(Jj 
{uj-ak)-Z R (k,Lu) 

|As(fc)| 2 ,„ s . 



(co + e(fc) + s fl (fc,-^)) 



\G b (k,u) = l (16) 



The gap function is assumed to have a d-wave form and 
it is related to the anomalous Green's function 

A s (fe) = J d 2 k'dujg(k - k')F(k',u) (17) 
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where g(k — k') is the d-wave pairing interaction. 

The strength of the superconductivity is determined 
by the magnitude of the gap squared. We assume a sim- 
ple parabolic form for T c (x)2& to mimic experiment (see 
FigE} and the gap is scaled with the superconducting, 
T c . 

We begin by examining the form of the LSR in 
the superconducting state. By inspecting the function 
G^ oh (k,0) defined in Eq.^Hl we see at once that the de- 
nominator diverges on the contours defined by £ (k) = 
and by £(fe) -I- £fl(fc,0) = 0. These are the same LSR 
contours that occur in the normal state. The only differ- 
ence is that the second set of contours now define zeroes 
of G^ oh (k, 0) not infinities. The Fermi surface of the hole 
pockets is gapped except along the nodal line where both 
A fl (fc) = A s (k) = 0. The form JTUJ for Gf o/l (fc, ui) then 
continues to satisfy the LSR. 

In the superconducting state the quasiparticle poles 
are given by solutions to a quartic equation 

(„ 2 - £ 2 (fc))(-^ 2 + Co(fe)) + 2A^(fc)(c 2 - umm 

+A|(fe)(o; 2 -^(fe))-A|(fe) = (18) 

which results in a further splitting of the quasiparticle 
bands. These are illustrated in Fig^| for a number of 
hole densities. The spectral weight redistribution is small 
when the original quasiparticle energies are away from 
the chemical potential. This can be seen for example in 
panel (d) of Fig 1 131 which shows the quasiparticle bands 
and their weight along the umklapp surface connecting 
(-7T, 0) — (0, 7r) . Here we see comparable weights only when 
the hole Fermi pockets are nearby. 

The superconducting energy gap modifies the DOS, 
Nt(ui). As illustrated in Fig^Jthe opening of the su- 
perconducting gap along the Fermi pockets (shown in 
Fig ll2fl leads to a pseudogap in Nt(uj) accompanied by 
van Hove singularities at the gap maxima. Note these van 
Hove singularities do not occur at the usual antinodal k- 



FIG. 12: (Color online) The dispersion (panels (a-c)) and 
weight (panel (d)) on the Luttinger surface with nonzero su- 
perconducting gap. There are 2 bands below the Fermi level, 
and shown here is the one closer to Fermi level. Around (n, 0) 
and (0, 7r), there is a substantial part of the spectral weight 
located at the lower band, only a small part remains on the 
band closer to the Fermi level as indicated in Fi dl3l 

points (±7r, 0), (0, ±7r) but rather on the extremities of 
the hole Fermi pockets. This is because the form Ea.|lUj 
we chose for the superconducting state does not repre- 
sent a merging of the two gaps A#(fc) and Ag(fc), rather 
both gaps keep their own identity. Further consideration 
of this point would require a better microscopic under- 
standing of the coexistence of d-wave superconductivity 
with the RVB spin liquid correlations than our simple 
phcnomcnological ansatz. 

Another clear feature in Nt(oj) is the van Hove singu- 
larity associated with antinodal saddle points in the RVB 
quasiparticle bands at negative energy. This feature does 
not appear in the STM experiments. This could be due 
to the strong local variations in the hole density. 



V. DISCUSSION AND CONCLUSION 

The ansatz for the Green's function Eq. [7j is motivated 
in the first place by the form of the Green's function de- 
rived recently by Konik and coworkersSiA^ for a doped 
spin liquid composed of an array of 2-leg Hubbard lad- 
ders. A second important input is the renormalizcd mean 
field derived many years ago by Zhang and coworkers 6 . 
As emphasized recently by Andseron and coworkers^, the 
RMF, although it treats the strong correlations simply by 
Gutzwiller renormalization factors, agrees qualitatively 
and even at times quantitatively with the results of varia- 
tional Monte Carlo (VMC) calculations using Gutzwiller 
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(a): (b): (c): (d): 

(0,0)-(ji,0) (n,0)-(n,n) (jt,ii)-(0,0) (7t,0)-(0,n) 

x=0.10 I I 1 I 

2 




FIG. 13: (Color online) Quasiparticle dispersion E(k) in the 
superconducting state CEq. llUI I along high symmetry direc- 
tions. The thickness of the lines is proportional to the spectral 
weight z(k)jg t . 



projected wavefunctions . The RMF and VMC calcu- 
lations at finite doping rely on a broken symmetry to 
introduce a gap at the Fermi energy. However, in the un- 
derdoped region the gap at the antinodal regions (~ J) 
is much larger than the expected superconducting tran- 
sition temperature T c as estimated for example by Wen 
and Leeii So one should expect this gap to persist in 
the normal state at temperatures T > T c , and so should 
be a property of the RVB spin liquid and not related 
to a broken symmetry. Another key feature of the un- 
derdoped region is that the Drude weight scales as the 
hole concentration. As a result the simplest explanation 
for these two properties is that the Fermi surface is par- 
tially truncated in a doped RVB spin liquid rather like 
the partial truncation that occurs through spin density 
waves in Cr and its alloys before the commensurate anti- 
ferromagnetic state is reached^. Here the key difference 
to the Cr alloys is the absence of a broken symmetry and 
long range order. Theses properties are reconciled by our 
ansatz for the normal state. In the superconducting state 
we introduced As and A# as separate gaps whereas in 
the RMF and VMC calculations there is only a single 




FIG. 14: (Color online) Density of states Nt(uj) in the pres- 
ence of a nonzero superconducting gap. The additional peaks 
closest to Fermi level are determined by the superconduct- 
ing gap which opens on the hole pockets for x < x c . For 
x > x c the peaks come from the antinodal directions for 

(±7T,0), (0,±7T). 



gap function opening up along a single Fermi surface. So 
it would appear that our ansatz is qualitatively different 
to the RMF and VMC theories in the superconducing 
state. However, one cannot be sure of this since the key 
quantity is the zero frequency Green's function G(fe,0) 
which is not directly available from a variational wave- 
function. We note that there have been recent reports 
on the numerical study of the pseudogap phase by us- 
ing cluster perturbation theory (CPT)2S and by the dy- 
namical mean field theory (DMFT)^£ of the 2-Dimension 
Hubbard model. 

Our ansatz does not deal with the origin of the d-wave 
superconductivity per sr. The relevant issue is the sta- 
bility of the normal state as a doped RVB state. This 
issue was addressed in the case of the doped array of 2- 
leg Hubbard ladders by KRT— . They showed that strong 
residual interactions acting on the hole pockets lead to a 
d-wave superconducting state or possibly a spin density 
wave state depending on the parameters. There is strong 
reason to believe that similar effects should occur here in 
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the 2-dimensional doped RVB spin liquid. Indeed such 
effect were conjectured by Honerkamp 8 , Laeuchli and 
coworkers^ in their analysis of the functional RG calcu- 
lation for the 2-dimensional t — t' — U Hubbard model. On 
a qualitative level we know that the RVB spin liquid has 
enhanced response in both d-wave pairing and antiferro- 
magnetic channels so that the Fermi pockets are moving 
in a background which is highly polarizable in both these 
channels which can lead to corresponding instabilities. 
The challenge is to develop a proper microscopic theory 
to describe the competition between the two instabilities. 
Note unlike the case of heavy fermions where a supercon- 



ducting dome straddles the QCP associated with onset 
of long range AF order—, here superconductivity seems 
to be most stable near the onset of the RVB phase with 
short range correlations. The good agreement that this 
ansatz displays with many properties of the anomalous 
normal state in the pseudogap region is evidence that it 
contains at least elements of the correct physics. 
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